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Abstract 

Epigenetic regulation plays essential role in cell differentiation and dedifferentiation, which are the in- 
trinsic processes involved in regeneration. To investigate the epigenetic basis of regeneration capacity, 
we choose DNA methylation as one of the most important epigenetic mechanisms and the MRL/MpJ 
mouse as a model of mammalian regeneration known to exhibit enhanced regeneration response in 
different organs. We report the comparative analysis of genomic DNA methylation profiles of the MRL/ 
MpJ and the control C57BL/6J mouse. Methylated DNA immunoprecipitation followed by microarray 
analysis using the Nimblegen '3 x 720 KCpG Island Plus RefSeq Promoter' platform was applied in order 
to carry out genome-wide DNA methylation profiling covering 20 404 promoter regions. We identified hun- 
dreds of hypo- and hypermethylated genes and CpG islands in the heart, liver, and spleen, and 3 7 of them 
in the three tissues. Decreased inter-tissue diversification and the shift of DNA methylation balance up- 
stream the genes distinguish the genomic methylation patterns of the MRL/MpJ mouse from the C5 7BL/ 
6J. Homeobox genes and a number of other genes involved in embryonic morphogenesis are significantly 
overrepresented among the genes hypomethylated in the MRL/MpJ mouse. These findings indicate 
that epigenetic patterning might be a likely molecular basis of regeneration capability in the MRL/MpJ 
mouse. 

Key words: MRL/MpJ mouse; regeneration; DNA methylation; epigenetics; genome-wide DNA methylation 
profiling 



1. Introduction 

The phenomenon of regeneration has become the 
object of extensive research in recent years. The investi- 
gations on the molecular basis of regeneration indiffer- 
ent models were rarely focused on epigenetic aspects, 
while it seems that these are the epigenetic mechan- 
isms which could play a key role in both maintaining 
the capacity for regeneration and remodelling involved 
in the regenerative response. All organs in an organism 



originate from a single cell (zygote), which clearly indi- 
cates that a single genome contains the information 
needed to grow organs and tissues, and possibly to 
restore lost and injured structures. This information is 
stored through the whole life. However, the expression 
of this information is gradually repressed during devel- 
opment, with epigenetic mechanisms including DNA 
methylation playing an essential part. 

Although the direct evidence on the involvement 
of DNA methylation in regeneration is scarce, a few 
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remarkable findings are worth mentioning including 
the role of the Shh gene methylation in amphibian 
limb regeneration 1 and the changes of DNA methyla- 
tion status in the Gadd45a promoter found during 
folate-induced regeneration response in the rodent 
spinal cord. 2 To our knowledge, no analysis of genome- 
wide DNA methylation profiles related to regeneration 
has been reported so far. 

The MRL/MpJ mouse, an inbred laboratory strain of 
Mus musculus, emerged as the model of mammalian re- 
generation in the 1 990s when the unusual phenom- 
enon of regenerative ear-hole closure was discovered 3 
to be followed by the observations of enhanced regen- 
erative response in other tissues including cornea, 4 
retina, 5,6 digit tips, 7,8 heart, 9 articular cartilage, 1 0 and 
spinal cord. The ability to close 2 mm through-and- 
through ear holes in the ear pinnae, which are used 
for a lifelong animal marking, clearly distinguishes the 
MRL/MpJ from most other mouse strains, although it 
isshared with the LG/J mouse,theMRL/MpJ's majoran- 
cestral strain. 1 1 The response to injury in the MRL/MpJ 
mouse does not result in scarformation, but in complete 
restoration of excised fragment with normal tissue archi- 
tecture, including blood vessels, hair follicles, sebaceous 
glands, and cartilage, that occurs within 30 days. 3 Since 
this process requires reinstatement of different types of 
tissues, it resembles epimorphic regeneration observed 
in amphibians. One important feature of injury response 
in the MRL/MpJ mouse is increased activity of matrix 
metalloproteinases that decompose the extracellular 
matrix breaking the basal lamina, a layer that initially 
covers the wound surface. This process allows the for- 
mation of the blastema-like structure of undifferenti- 
ated cells that drive the process of regeneration. 12 
Another important feature of the MRL/MpJ mouse is 
its abnormal cell cycle profile with accumulation of 
the cells in the G2/M phase. This hallmark has been 
connected with the deficiency of p2 1 , an inhibitor of 
cell cycle progression, that is not detected in the 
tissues of the MRL/MpJ mouse even after 7-irradiation. 
Although the regeneration abilities of the MRL/MpJ 
expand themselves on other tissues, they do not reach 
as spectacular effects as those observed in urodele 
amphibians which restore lost limbs 13 and regenerate 
as serious injuries as the complete transection of the 
spinal cord 14 and partial heart resection. 1 5 In the MRL/ 
MpJ mouse, lost digits do not re-grow, 8 heart regener- 
ation is rather limited, 9,16,17 and the heart does not 
heal in some types of injuries. 1 8-22 The regenerative 
phenotype of the MRL/MpJ mouse has been found to 
be multigenic. 23 It is also remarkable that increased re- 
generative response in the MRL/MpJ mouse is found in 
different tissues, and the strain maintains selected fea- 
tures of embryonic metabolism in adults including the 
enhanced expression of key stem cell markers, Sox2 and 
Nanog 24 ' 25 



The object of the presented study is the comparative 
analysisof the genome-wide DNA methylation profiles 
in the MRL/MpJ "healer" mouse, using the C57BL/6J 
mouse as the reference. In our attempt to understand 
the epigenetic basis of regeneration, we would like 
to identify the differentially methylated genes and to 
evaluate their possible association with the regenera- 
tive capacity of the MRL/MpJ mouse. The study was 
designed to analyse global DNA methylation profiles 
in heart tissue, as it is known to have poor, if any, regen- 
erative capacity, in order to refer to the liver, which 
has a strong regenerative potential, in addition to 
spleen, selected for its involvement in cell-mediated 
pathways of the immune system. As the regenerating 
tissues are difficult to obtain in sufficient amounts 
and the DNA methylation patterns in cultured cell 
lines are prone to change, we would like to focus 
the study on the normal tissues, thus examining the 
basis of regeneration rather than the regenerative 
response. 



2. Materials and methods 

2.1 . Tissue samples and nucleic acid extraction 
Tissue samples (hearts, livers, and spleens) from 

5-week-old males and 8-week-old females of MRL/MpJ 
(stock #000486) and C57BL/6J (stock #000664) 
murine strains were purchased from the Jackson 
Laboratory. The samples from 8 week-old females 
were not used for the microarray analyses, but for 
the validation experiments only. Genomic DNA was 
isolated using the DNeasy Blood &Tissue Kit (QIAGEN, 
cat. no. 69504) and RNA with the RNeasy Mini Kit 
(QIAGEN, cat. no. 741 04) according to the manufac- 
turer's protocols. Tissues were delivered in RNAlater 
stablilization reagent (QIAGEN, cat. no. 76104) and 
were disrupted in liquid nitrogen before nucleic acids 
isolation. 

2.2. Methylated DNA immunoprecipitation 
Methylated DNA immunoprecipitation (MeDIP) was 

performed according to the modified Weber's proto- 
col. 26 Equal amounts of genomic DNA from three indi- 
viduals were pooled for each strain and tissue. Each 
pooled sample (7 |xg) was diluted in 1 50 |xl of water 
and randomly sheared by sonication to generate frag- 
ments between 300 and 1000 bp. A 4-|xg portion 
of sonicated DNA was diluted in 450 julL of 1 0 mM 
Tris-HCl, pH 8.0, heat denatured at 98°C for 1 0 min 
and immediately cooled on ice for 1 0 min. A portion 
of the untreated, sonicated DNA was left to serve 
as the input control. The CpG-methylated DNA frag- 
ments were precipitated with anti-5-methyl cytidine 
antibody (Eurogentec cat. no. BI-MECY-1 000), and 
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the antibody-DNA complexes were captured with 
Dynabeads (Invitrogen, cat. no. 1 1 2 02). The collected 
beads were washed in order to remove non-specifical- 
ly-bound DNA,followingthetreatmentwith proteinase 
Kfor 20-24 h at 50°C in orderto remove proteins. The 
CpG-methylated DNA was extracted with the phenol:- 
chloroform method, precipitated with ethanol and 
glycogen, and resuspended in 60 jxl of 1 0 mM Tris- 
Cl, pH8.0. 

2.3. DNA Labelling and Hybridization 

DNA labelling and hybridization was performed 
according to NimbleGen's protocol with our modifica- 
tions. 27 The immunoprecipitated CpG-methylated 
DNA (test) and the untreated, sonicated DNA (input 
control) were labelled by using the random priming 
with the Nimblegen Dual-Color DNA Labelling kit 
(Roche-Nimblegen, cat. no. 06370250001) with 
fluorescent dyes Cy3 and Cy5, respectively. The com- 
bined test and input DNA samples were suspended in 
hybridization buffer (Roche, cat. no. 05583683001) 
co-hybridized onto Mouse DNA Methylation 3 x 
72 0 K CpG Island Plus RefSeq Promoter Arrays for 
20 h at 42°C, following washing with the Nimblegen 
Wash Kit (Roche cat. no. 05584507001 ). 

2.4. Microarray data acquisition and processing 

We performed image acquisition with an MS200 
Scanner (Roche, NimbleGen) at 2 |jim resolution by 
using high-sensitivity and autogain settings. The 
data from scanned images were extracted and pro- 
cessed with DEVA v. 1.0.2 (Roche, Nimblegen) using 
default parameters. Data processing included obtain- 
ing log 2 ratios, P-scores, and peak identification. Log 2 
ratios represent the ratios of the immunoprecipitated 
DNA signal to the input DNA signal. The mean log 2 
ratios were calculated using a Tukey biweight mean. 
P-scores were calculated from log 2 ratios by performing 
the sliding-window (7 50 bp) Kolmogorov-Smirnov 
(KS) test around each probe applied to determine 
a -log1 0 P-value from the windowed KS test around 
that probe. The methylation peaks were obtained 
by merging consecutive probes with P-scores over 
the assumed cut-off (default 2.0), with maximum 
500 bp spacing and minimum of two probes per 
peak. The processed data (P-scores) and raw data 
were deposited in the Gene Expression Omnibus 
(GEO) database under the accession number 
GSE49221. 

2.5. Identification of differentially methylated genes 
The methylation peaks were mapped to features 

(transcription start sites, primary transcripts, CpG 
islands, and other tiled regions) using DEVA v. 1 .0.2, as- 
suming the default distances 5000 bp upstream and 



1 000 bp downstream the feature The annotated data 
for the methylation peaks obtained at cut-off of 1 .0 
and 2.0 were merged using GenixNet, a script created 
by Dr P. Madanecki and adjusted to the analysed data- 
sets, so that a single table containing all genes, tiled 
regions, and CpG islands with the methylation peak 
values for each strain and tissue at two different cut- 
off points was created. The data merged using 
GenixNet were sorted to identify the differentially 
methylated genes using the Excel spreadsheet. If a 
peak was mapped to more than one feature (e.g. tran- 
scription start site and primary transcript), the data 
listed in the tables of differentially methylated peaks 
(Supplementary Table S1) were restricted to a single 
feature selected according to the following order: 
transcription start site, primary transcript, CpG island, 
and tiled region. 

2.6. Validation of the MeDIP/ microarray results 
and gene expression level quantification 

The methylation status of selected CpG nucleotides 
was evaluated using bisulphite sequencing and/or 
genomic DNA digestion with a CpG methylation- 
sensitive HpaW restriction endonuclease. Gene expres- 
sion levels were quantified using real-time PCR. 

The bisulphite conversion of genomic DNA was 
carried out with the EZ DNA Methylation™ Kit (Zymo 
Research, cat. no. D5002) using 500 ng of genomic 
DNA. Approximately 1 00 ng of obtained converted 
DNA was used for the subsequent PCR reaction with 
the Maxima Hot Start Taq DNA polymerase system 
(ThermoScientificBio, cat. no. EP0601). The PCR reac- 
tions were performed on Veriti® 96-well Thermal 
Cycler (Applied Biosystems). The PCR amplified DNA 
fragments were extracted from the agarose gel and 
purified using the ISOLATE Gel and PCR kit (Bioline, 
cat. no. BIO-52029).The PCR products were sequenced 
using the Applied Biosystems ABI 3730XL/ABI3700 
DNA sequencer (Genomed, Poland), and the chroma- 
tograms were analysed with ContigExpress (VectorNTI 
1 1 .0, Invitrogen). 

The CpG methylation-sensitive digestion of genomic 
DNA was carried out with the EpiJET DNA Methylation 
Analysis Kit based on Msp\/Hpa\\ digestion (Thermo- 
ScientificBio, cat. no. K1441) using 500 ng of genomic 
DNA and 50 ng of digested DNA for quantitative real- 
time PCR analysis. The methylation level was quantified 
using the 2-ACt method and expressed as the ratio 
between /-/pfl//-digested DNA (target) and input/non- 
digested DNA (reference) used for enzymatic reaction. 
cDNA synthesis was performed with Maxima Reverse 
Transcriptase (ThermoScientificBio, cat. no. EP0742) 
using 200 ng of total RNA and oligodT 20 . Approximately 
1 0 ng of cDNA was used for subsequent real-time PCR 
reactions. The gene expression levels were calculated 



608 



Genomic DNA Methylation Profiles in the MRL/MpJ 



[Vol. 20, 



usingthe 2-ACt method \N\tbActb as the reference gene. 
The results are presented as means + SD. 

Real-time PCR reactions were carried out with 
FastStart Essential DNA Green Master (Roche, cat. 
no. 06402712001) on LightCycler® Nano (Roche). 
The primer sequences are listed in Supplementary 
Table S2. 

2.7 . Linear regression and scatter plots 

Linear regression coefficients were determined, and 
scatter plots were made using the XLSTAT package 
(Addinsoft). 

2.8. Chromosomal ideograms with the locations of 
differentially methylated genes 

The chromosomal mapswiththe locations of differen- 
tially methylated genes were plotted using Idiographica 
v. 2.1 (http://www.ncrna.org/idiographica). 28 

2.9. Counting the DNA methylation peaks at different 
value ranges 

The lists of unique peaks were obtained from the lists 
of mapped peaks generated by DEVA v. 1 .0.2. The 
numbers of peaks at different methylation values 
were calculated in Excel spreadsheet. 

2. 1 0. Gene set enrichment analysis 

The gene set enrichment analyses were performed 
and functional annotations were assigned by using 
DAVID Bioinformatics Resource 6.7 [National Institute 
of Allergy and Infectious Diseases (NIAID), National 
Institute of Health 29 ' 30 ]. An EASE Score, a modified 
Fisher's Exact test, was applied to calculate statistical sig- 
nificance of gene set enrichments. 

2.11. Processing and analysis of genomic DNA 
methylation microarray data downloaded from 
the Gene Omnibus Expression Database 

2.1 1.1. Gene expression data The microarray data 
on gene expression in the heart and liver of the MRL/ 
MpJ mouse were retrieved from GSE1 9322,GSE471 0, 
GSE25322, and GSE5241 datasets deposited in the 
GEO database. The gene expression values were 
merged with the annotations from the corresponding 
platform files in Excel spreadsheet. The logarithmed 
signal values were exponentiated. Average values were 
calculated from replicates. To calculate median and 
quartile expression values for selected gene subsets 
(e.g. the genes hypomethylated in the MRL/MpJ 
mouse), the expression values were extracted in Excel 
spreadsheet using gene identifiers. Gene expression 
va I ues between the su bsets of genes se lected accord i ng 
totheir methylation levels were compared by Wilcoxon 
rank-sum two-tailed tests using the SAS statistical 
package version 9.3 (SAS Institute, Inc., Cary, NC, USA). 



2.11.2. Gene methylation data Processed DNA 
methylation data (P-scores) were downloaded from 
theGSE21 41 5 dataset deposited in the GEO database. 
The genome coordinates were converted from the NCBI 
MM8 build to the NCBI MM9 build using 'Lift Genome 
Annotations' service from the University of California 
Santa Cruz Genome browser (http://genome.ucsc. 
edu/cgi-bin/hgLiftOver). The methylation peaks were 
obtained with NimbleScan 2.6 (Roche, Nimblegen) at 
cut-off of 1.0 using default parameters (at least two 
probes per peak with maximum spacing of 500 bp). 
Average values were calculated from replicates. The 
methylation peaks were mapped to gene features 
using DEVA v. 1.0.2, assuming the standard distance 
5000 bp upstream and 1 000 bp downstream the 
feature. Average values were used for the features with 
more than one methylation peak. 

2.12. Analysis of correlations between DNA methylation 
and gene expression data 

Gene expression and methylation data were merged 
using GenixNet. Pearson's correlation coefficients 
were calculated using the SAS statistical package 
version 9.3 (SAS Institute, Inc.). 

2.13. Analysis of correlations between genome-wide 
DNA methylation profiles 

As the formats of microarray platforms applied in this 
study and that by Liang etal. are different, we restricted 
the analysis to the peaks mapped to the transcription 
start sites shared by both platforms (1 8 1 99 unique 
genes included in Roche MM8 3 85 K RefSeq Promoter 
Array). If more than a single peak was mapped to 
a gene, average values were calculated by using 
Excel 'pivot tables' function. The methylation peaks 
for this analysis were obtained at cut-off of 1.0 in 
order to reduce the number of '0' under the cut-off, 
which, in fact, represent a spectrum of values. The 
peaks under the cut-off of 1 .0 were given the value '0'. 
The arrays were merged using Entrez identifiers. 
Pearson's correlation coefficients were calculated 
using the SAS statistical package version 9.3 (SAS 
Institute, Inc.). 



3. Results and discussion 

We investigated the epigenetic basis of the regener- 
ation phenomenon in the MRL/MpJ mouse using 
genome-wide DNA methylation profiling. As the MRL/ 
MpJ strain does not have an immediate relative, which 
does not exhibit enhanced regeneration abilities, we 
chose the C57BL/6J mouse as the reference strain. We 
carried out genome-wide DNA methylation profiling 
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in hearts, livers, and spleens of the MRL/MpJ and 
C57BL/6J strains using the MeDIP approach followed 
by microarray analysis. The Nimblegen microarray 
system '3 x 720 K CpG Island Plus RefSeq Promoter' 
we applied in the study includes 20 404 promoter 
regions; 22 881 transcripts and 1 5 980 CpG Islands. 
The DNA methylation profiles we obtained for the 
MRL/MpJ and C57BL/6J were compared in order to 
identify the differentially methylated genes and 
regions. The overall comparisons of inter-strain similar- 
ities and differences between the examined genomic 
DNA methylation profiles are visualized using scatter 
plots (Supplementary Fig. S1 ). 

This study was designed as an initial analysis of 
genome-wide DNA methylation patterns and we 
restricted the selection of tissues to heart, liver, and 
spleen. Heart and liver were selected to contrast the 
non-regenerating organ with the regenerating one. 
Spleen was chosen for its role in the immune 
system, as immune response to injury seems to be 
one of the most critical events in the regeneration 
process. 31 



3.1 . Genes differentially methylated in the MRL /MpJ 
mouse 

As the DNA methylation patterns are prone to indi- 
vidual and conditional variation, we focused on the ana- 
lysis of the genes revealing the most distinguishing 
differences in the methylation status. We assumed a 
gene/region to be differentially methylated if a methy- 
lation peak mapped to this gene/region meets at least 
one of two conditions: either the peak value is at least 
3.0 or higher and it is at least twice higher than in the 
reference (at cut-off of 1 .0) or it is at least 2.0 and it is 
lower than 1 .0 for the reference (see Section 2 for the 
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Figure 1. The numbers of differentially methylated genes and CpG 
islands in the MRL/MpJ mouse in the heart, liver, and spleen 
tissues. The numbers of CpGs distant from the genes (out of 
- 5000 bp and + 1 000 bp range) in brackets. It is assumed here 
that a gene/region is differentially methylated if it meets at least 
one of two conditions: either (i) the peak value is at least 3.0 or 
higher and it is at least twice higher than in the reference or (ii) 
it is at least 2.0 and it is lower than 1 .0 for the reference. 



definition of peak value). A peak is mapped to a 
feature (transcription start site, primary transcript, 
CpG island, or tiled region) if it is located within 
5000 bp upstream and 1000 bp downstream the 
region (a transcription start site, primary transcript, or 
CpG island). 

As the result several hundreds of differentially 
methylated genes, we will refer to as 'hypomethylated' 
and 'hypermethylated', were identified in liver, spleen, 
and heart tissue of the MRL/MpJ mouse (Fig. 1 and 
Supplementary TableS! ).The majority of differentially 
methylated genes and regions were not shared by 
the three examined tissues: heart, liver, and spleen. 
Twenty-three genes and CpG islands were found to be 
hypomethylated and 1 4 to be hypermethylated in 
the MRL/MpJ mouse in the three examined tissues. 
These genes and CpG islands were located on all 
chromosomes except the 1 8, 1 9, X, and Y (Fig. 2), and 
they are associated with a number of different 
functions. The complete list and the descriptions 
of these genes are presented in Supplementary 
Table S1 , while a selection of remarkable representa- 
tives of this group (Table 1) are discussed later in this 
article. 

3.2. Validation of the MeDIP /microarray results 

A selection of DNA methylation peaks were validated 
in the MRL/MpJ and the C57BL/6J reference using 
bisulphite sequencing and/or digestion with a CpG 
methylation-sensitive restriction endonuclease HpaW, 
followed by quantitative real-time PCR analysis [Supple- 
mentary Table S3 (summarized results), Supplementary 
File S1 (experiment details), and Supplementary File S2 
(representative DNA sequencing chromatograms 
assembled with the reference sequences)]. Gene expres- 
sion levels of selected genes were examined by quantita- 
tive real-time PCR analysis [Supplementary Table S3 
(summarized results) and Supplementary File S1 (experi- 
ment details)]. 

3.3. The chromosomal locations of differentially 
methylated regions 

The hyper- and hypomethylated genes and regions 
were found on all chromosomes. Several chromosomal 
loci display significant positional enrichments of the 
genes differentially methylated in the MRL/MpJ 
mouse, but no remarkable regions of peak accumula- 
tion could be distinguished (Fig. 2A and B). 

3.4. Gene methylation status and gene expression 
DNA methylation of gene promoters is known to 

repress gene expression. With regard to evaluating the 
impact of DNA methylation on gene expression in the 
examined tissues, we analysed the available expression 
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Figure 2. The chromosomal locations of genes and CpG islands differentially methylated in the MRL/MpJ mouse. (A) The genes and CpG islands 
that are differentially methylated in the three examined tissues: heart, liver, and spleen. PT- DNA methylation peak mapped to a primary 
transcript, kb — distance in kb is given if the methylation peak is located out of the range 5000 bp upstream and 1 000 bp downstream 
the feature. (B) The genes and CpG islands differentially methylated in the heart, liver, and spleen. The genes displaying substantially 
higher and lower methylation degree in the MRL/MpJ mouse are marked red (hypermethylated, grey in the print) and blue 
(hypomethylated, black in the print), respectively. The values and genomic coordinates of the DNA methylation peaks as well as the 
descriptions of mapped of genes are listed in Supplementary Table S1 . Heart vs. liver analysis singled out the genes and CpG islands, 
which have different methylation status in the heart than in the liver of the MRL/MpJ, but not the C57BL/6J, mouse, (see the chapter: 
'Genes with similar DNA methylation status in the heart and liver of the MRL/MpJ but not the C57BL/6J mouse')- 



Table 1. The remarkable representatives of the genes differentially methylated in the MRL/MpJ mouse P 
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Gene 


Full name 


Peak location 


Methylation 
peak MRL/MpJ 


Methylation 
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Phc3 


Polyhomeolic-like 3 (Drosophila) 


TSS 


2.3 


2.8 


2.1 


0.0 


0.0 


0.0 


Polycomb group (PcG) ; transcriptional repression of many 
genes including Hox genes 


Fggy 


FGGY carbohydrate kinase domain containing 
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4.4 
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Carbohydrate kinase (E) 


Gcnlll 


GCN1 general control of arnfno-acid synthesis 1-like 1 (yeast) 


TSS 


2.8 


3.4 


3.2 


0.0 


1.5 


1.4 


A possible activator of Ras super-family of small GTPases 
(E) 


Rnf7 IRasal 


Ring finger protein 7/RAS p2 1 protein activator 2 


58 kb from TSS 
12 kb from 3' 




4.4 


4.1 


4.0 


1.0 


1.4 


1.5 


Proteasomal degradation; cell cycle progression/inhibitory 
regulator of the Ras-cAMP pathway 


Tbrgl 


Transforming growth factor beta- regulated gene 1 


TSS 


2.5 


2.2 


2.7 


0.0 


0.0 


0.0 


Growth inhibitor; TP53/p53 activation, causes Gl arrest 


Etfa 


Electron transferring flavoprotein. alpha polypeptide 


TSS 


2.6 


2.4 


2.1 


0.0 


0.0 


0.0 


Electron transfer flavoprotein 


PpiM 


Peptidylprolyl isomerase (cyclophilin)-like 4 


TSS 


2.7 


3.3 


3.3 


0.0 


1.2 


0.0 


Protein folding 


Nlrplb 


NLR family, pyrin domain containing IB 


TSS 


2.8 


2.1 


3.0 


0.0 


0.0 


0.0 


Mediator of macrophage susceptibility 


Nlrplc-ps 


NLR family, pyrin domain containing 1C. pseudogene 


TSS 


2.7 


2.3 


2.4 


0.0 


0.0 


0.0 


Inflammation mediator 


Grebl 


Gene regulated by oestrogen in breast cancer protein 


TSS 


2.1 


2.9 


2.4 


0.0 


0.0 


0.0 


Critically involved in the oestrogen-induced growth of breast 
cancer cells 32 


Akrlel 


Aldo-keto reductase family 1, memberEl 


TSS 


4.6 


3.5 


3.6 


0.0 


0.0 


0.0 


Synthesis of 1,5-anhydro-D-glucitol. 


Tob2 


Transducer of ERBB2, 2 


TSS 


4.1 


3.0 


3.6 


0.0 


0.0 


0.0 


Inhibits cell cycle progression from the G0/G1 to S phases 


Pdcd2 /Prdm9 


Programmed cell death 2/PR domain containing 9 


9 kb from TSS 
/6 kb from 3' 


4.3 


2.3 


4.7 


0.0 


0.0 


0.0 


Haematopoietic stem cell differentiation 1 Vspecific 

methylation of 'Lys-4' of histone H3 


2300002M23Rik 


RIKEN cDNA 2300002M23 gene 


TSS 


2.8 


3.5 


3.2 


0.0 


1.3 


1.1 


Interacts with extracellular matrix proteins 


Tfg 


Trk-fused gene 


TSS 


1.0 


0.0 


1.2 


3.1 


3.9 


3.5 


Fusion oncoprotein 


Atp6vla 


ATPase. H+ transporting, lysosomal VI subunit A 


TSS 


0.0 


0.0 


1.7 


3.1 


2.9 


3.8 


Acidifying of intracellular compartments 


NaaSO 


rV(alpha)-acelyllransferase 50. NatE catalytic subunit 


TSS 


0.0 


0.0 


1.7 


3.1 


2.9 


3.8 


Protein acetylation; chromosome segregation. 


KirrelS 


KinofIRRElike3 


760 kb from TSS 


0.0 


0.0 


1.7 


3.0 


2.6 


3.8 


Responsible for supporting hematopoietic stem cells 


Ybx2 


Y-box protein 2 


TSS 


0.0 


0.0 


1.2 


2.7 


3.1 


3.9 


Regulation of the stability and/or translation of germ cell 
mRNAs 


Fa2h 


Fatty acid Ihydroxylase 


TSS 


0.0 


0.0 


1.4 


2.3 


4.2 


3.6 


Formation of alpha-hydroxylated sphingolipids 


Ntng2 


Netrin G2 


0.14 kb from 3' 


0.0 


1.2 




0.0 


2.4 


3.1 


3.5 


Promotes neurite outgrowth 


CnotS 


CCR4-NOT transcription complex, subunit 8 


TSS 


0.0 


0.0 


0.0 


2.9 


3.7 


3.2 


Ubiquitous transcription factor 


Cacnala 


calcium channel, voltage-dependent. P/Q type, alpha 1 A subunit 


PT 


0.0 


0.0 


0.0 


2.4 


4.0 


2.7 


Calcium channel associated with migraine in humans 


Asapl 


ArfGAP with SH# domain, ankyrin repeal and PH domainl 


PT 


0.0 


0.0 


0.0 


2.4 


3.5 


2.9 


May function as a signal transduction protein involved in 
the differentiation of fibroblasts into adipocytes 


TptI 


Tumour protein, translationally controlled 1 


TSS 


0.0 


0.0 


0.0 


2.2 


3.4 


3.2 


Calcium binding and microtubule stabilization 


Chrdl2 


Chordinlike 2 


26 kb from TSS 


0.0 


0.0 


0.0 


2.9 


2.7 


2.7 


Negative regulator of cartilage formation; interacts with 
bone morphogenetic proteins 


Foxc2 


Forkhead box C2 


PT 


0.0 


0.0 


0.0 


2.2 


3.3 


2.5 


Transcriptional activator expressed temporally in developing 
embryos 


Insig2 


Insulin-induced gene 2 


TSS 


0.0 


0.0 


0.0 


2.7 


2.4 


2.7 


Control of cholesterol synthesis 


Fbxo41 


F-box protein 41 


PT 


0.0 


0.0 


0.0 


2.6 


2.8 


2.3 


E3 ubiquitin ligase complex 


SmydS 


SET and MYND domain containing 5 


TSS 


0.0 


0.0 


0.0 


2.2 


2.5 


2.9 


Up-regulated by retinoic acid treatment in embryonic 
carcinoma cells 


LrrcSd 


Leucine rich repeal containing 8D 


TSS 


0.0 


0.0 


0.0 


2.7 


2.3 


2.3 


Up-regulated after cell cycle exit 34 



<Tj' 



Prpf39 


PRP39 pre-mRNA processing factor 39 homologue ( yeast) 


TSS 


0.0 


0.0 


0.0 


2.2 


2.8 


2.3 


Pre-mRNA splicing 


Sk4a4 


Solute carrier family 4 (anion exchanger), member 4 


TSS 


0.0 


0.0 


0.0 


2.2 


3.4 


3.1 


Involved in the regulation of bicarbonate secretion and 
absorption and intracellular pH (E) 


Tc2n 


Tandem C2 domains, nuclear 


TSS 


0.0 


0.0 


0.0 


2.1 


2.7 


2.4 


Possible association with elevated levels of factor VIII and 
von Willebrand Factor 35 


Naa30 


N{ alpha )-aeetyllransferase 30, NatC catalytic subunit 


TSS 


0.0 


0.0 


0.0 


2.2 


2.7 


2.3 


Probable A'-acetyltransferase. 


Zcchc9 


Zinc finger. CCHC domain containing 9 


TSS 


0.0 


0.0 


0.0 


2.4 


2.1 


2.3 


Suppresses the MAPK signal pathway 


Jakinipl/Crmpl 


Janus kinase and microtubule interacting protein 1/collapsin response 
mediator protein 1 


38 kb from 3' /78 kb 
from TSS 


0.0 


1.4 


0.0 


3.0 


2.9 


3.2 


Microtubule-dependent transport of the GAB A-B 
receptor/remodelling of cyloskeleton 


Sp9 


Trans-acting transcription factor 9 


PT 


0.0 


0.0 


0.0 


2.4 


3.0 


2.4 


FgflO and Wnt/beta-catenin signaling during vertebrate 
limb outgrowth 


Mbd3 


Methyl-CpG binding domain protein 3 


TSS 


1.2 


0.0 


1.1 


2.6 


3.5 


3.2 


Recruits histone deacetylases and DNA melhyltransferases; 
gene silencing 


Shh 


Sonic hedgehog 


PT 


0.0 


0.0 


0.0 


1.2 


2.5 


2.4 


Embryonic patterning and morphogenesis 


Gatadla 


GATA zinc finger domain containing 2B 


TSS 


0.0 


0.0 


0.0 


2.2 


0.0 


0.0 


Transcriptional repressor; complexes with MBD2, HDAC1 , 
HDAC2 


Epha2 


Eph receptor A2 


9.3 kb from TSS 


0.0 


0.0 


0.0 


2.0 


0.0 


2.5 


Post-natal angiogenesis 


Esrrg 


Oestrogen-related receptor gamma 


390 kb from TSS 


0.0 


0.0 


0.0 


3.8 


0.0 


1.3 


Transcription activator 


DrmtSb 


DNA methyltransferase 3B 


TSS 


0.0 




1.3 


2.2 


2.5 


2.0 


De novo melhylation; development; expressed in ESC and 
ICM 


Pax! 


Paired box gene 2 


22 kb from TSS 


0.0 


0.0 


1.7 


2.0 


0.0 


2.4 


Embryo development 


Gadd4Sa 


Growth arrest and DNA damage-inducible 45 alpha 


TSS 


2.3 


0.0 


2.4 


0.0 


1.3 


0.0 


Stimulates DNA repair; inhibits cell entry into S phase 


Abliml 


Actin-binding LIM protein 1 


TSS 


2.3 


3.5 


4.1 


0.0 


2.8 


4.5 


LEVI domain; development 


Fhll 


Four and a half LIM domains 1 


TSS 


3.1 


2.6 


3.3 


1.5 


2.1 


2.2 


LIM domain: muscle development or hypertrophy 


Csrpl 


Cysteine and glycine-rich protein 1 


TSS 


2.9 


2.6 


2.4 


1.3 


2.3 


0.0 


LEVI domain; neuronal development 


Hdacl 


Histone deacetylase 1 


TSS 


2.9 


2.2 
0.0 


2.5 


1.1 


2.4 


1.5 


Epigenetic repression; cell cycle progressionand 
development 


Anxal 


Annexin Al 


TSS 


2.4 


1.4 


0.0 


0.0 


1.1 


Cardioprotecuon 


Gngtl 


Guanine nucleotide binding protein (G protein) 


TSS 


2.1 


1.2 


1.4 


0.0 


2.3 


0.0 


Cardioprotection 


Pik3rl 


Phosphatidviinositol 3-kinase. regulatory subunit, polypeptide 1 


TSS 


3.6 




1.7 


1.3 


0.0 


1.4 


Cardioprotecuon 


MapkS 


Mitogen-activated protein kinase 8 


TSS 


0.0 


0.0 


0.0 


0.0 


3.0 


0.0 


Cell signalling 


Map4k5 


Mitogen-aclivaled protein kinase kinase kinase kinase 5 


TSS 


1.2 


0.0 


1.1 


2.8 


4.0 


3.3 


Cell signalling 


Smad3/Smad6 


MAD homologue 3 (Drosophila)/MAD homologue 6 (Drosophila) 


170 kb from 3' 
/ 25 kb from TSS 


1.2 


1.2 


4.2 


0.0 


0.0 


1.7 


TGF-beta signalling 


Pdgfa 


Platelet-derived growth factor, alpha 


TSS 


0.0 


1.3 


0.0 


1.3 


2.3 


2.3 


Potent mitogen; released upon wounding stimulates 
adjacent cells to grow 


Pdgfra 


Platelet derived growth factor receptor, alpha polypeptide 


TSS 


0.0 


0.0 


1.3 


2.1 


1.5 


2.7 


Receptor of Pdgfb and Pdgfa; cell proliferation and 
development 


Pdgfrl 


Platelet-derived growth factor receptor-like 


TSS 


2.7 


0.0 


2.5 


2.2 


2.8 


2.3 


Homologue to the ligand binding domain of Pdgfrb, 
associated with Behcet disease 36 



Scale 


0.0 


1.0 


2.0 


3.0 


4.0 


5.0 1 

















If not indicated otherwise, the information on gene functions is taken from DAVID and Gene Entrez (E, http://www.ncbi.nlm.nih.gov/gene/). 
H: heart; L: liver; S: spleen; TSS: transcription start site; PT: primary transcript. 
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Gene Expression for Different Gene Methylation Ranges. 
Mouse Male Hearts 8 wks. 



a «. 

I 80 
| 60 
| 40 
20 



| 

a 
i 

5 
2 

MRL/Mpj dataset: GSE4710 



methylation 
peak value 
range 



CS7BL/6J dataset: GSE4710 



Gene Expression for Different Gene Methylation Ranges. 
Mouse Livers 15-20 wks, chow diet, fasted 4 h prior collection. 
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Gene Expression for Different Gene Methylation Ranges. 
Mouse Male Hearts 15-20 wks 
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Gene Expression for Different Gene Methylation Ranges. 
Mouse Livers, 11 wks, injected 0.1 ml corn oil twice a week for 42 days 
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GSE25322 



p-value 


0.0168 


0.0002 


0.0308 


0.0013 
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< 0.0001 


< 0.0001 



Figure 3. Gene methylation and gene expression profiles. (A) Gene expression levels at different DNA methylation peak value ranges in heart 
and livers of the MRL/MpJ mouse. The microarray gene expression data representing the tissues collected from uninjured organs were 
retrieved from the following datasets GSE471 0 and GSE1 9222 (hearts), GSE5241 and GSE25322 (livers) that had been obtained in 
four independent studies. 1 ^' 7,37,38 Wilcoxon rank-sum two-tailed tests were performed in order to compare gene expression values 
between the groups of genes selected according to the DNA methylation level ranges (determined in this study). The median and 
quartile values were normalized to the highest value from each dataset. (B) The correlations between the DNA methylation and gene 
expression levels. The Pearson's correlation coefficients were calculated for the gene methylation values determined in this study for the 
heart and liver of the MRL/MpJ and C57BL/6J mouse, and corresponding gene expression values that were obtained from four other 
independent studies, 1 6,1 7 - 37 > 38 as described above. 



microarraydatafortheMRL/MpJ heartand livertissues 
reported in four different studies. 16,1 7,37,38 Our evalu- 
ation was based on the correlations between DNA 
methylation status and expression signals, as well as 
mean values of expression signals calculated for 



subsets of genes at different peak value ranges. We 
found that the expression signal values were increasing 
with the decrease of DNA methylation status (Fig. 3A), 
and the values of methylation peaks and gene expres- 
sion signals were reversely correlated (Fig. 3B). The 
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Pearson's correlation coefficients for gene methylation 
and expression values in the heart though significant for 
both theC57BL/6J and MRL/MpJ mouse were lowerfor 
the latterone. This analysis confirmed that the obtained 
genomic DNA methylation patterns a re adequate to the 
transcriptomic profiles. However, the gene expression 
signals in hearts for the groups of genes differentially 
methylated in the MRL/MpJ and C57BL/6J did not 
show statistically significant differences between 
these two strains (Fig. 3A). This inconsistence could 
be explained by age differences, as we used DNA methy- 
lation data obtained for 5-week-old animals and the 
gene expression data had been collected for much 
older mice (8- to 20-week old). It should also be 
noted that the majority of DNA methylation peaks 
mapped to the differentially methylated genes in the 
MRL/MpJ are located at higher distance than those in 
the C57BL/6J mouse (Fig. 6B). After all, in spite of the 
age differences, the overall correlations between DNA 
methylation and gene expression were confirmed. 

We also would I ike to punctuate three possible reasons 
why gene methylation and expression levels could be in- 
consistent. First, it should be noted that tissues unlike cell 
lines are more likely to be mosaics of cells that show dif- 
ferent DNA methylation and gene expression patterns, 
thus a methylation peak may represent gene methyla- 
tion status in a fraction of tissue cells, while another frac- 
tion of cells is responsible for high gene expression. 
Secondly, an individual gene does not have to show 
higher expression when unmethylated, as other repres- 
sion mechanisms could be involved. Thirdly, some 
genes are induced due to demethylation of selected 
CpGs within a CpG island. 39 

3.5. Functional profiling of genes differentially 
methylated in the MRL/MpJ mouse: gene set 
enrichment analysis 
Once the genesdifferentially methylated in the MRL/ 
MpJ mouse were singled out, we could examine which 
functional and structural categories were overrepre- 
sented among these genes, to be followed by looking 
for associations with the phenotype. As the first step, 
we carried out gene set enrichment analysis for the 
groups of hypo- and hypermethylated genes using 
the DAVID algorithm and database. 29,30 We examined 
the sets of genes hypo- and hypermethylated in hearts, 
livers, and spleens (six gene sets, Supplementary 
Table S1). The analysis revealed a numberof statistically 
significant functional terms (Supplementary Table S4). 
Although terms which are not statistically significant, 
as well as those found for one of tissues, could provide 
useful clues on the association between DNA methyla- 
tion and regeneration potential, we decided to focus 
on the statistically significant terms representative of 
the three examined tissues. The statistically significant 



functional terms representing the hypomethylated 
genes in the three examined tissues (Supplementary 
Table S5) could be associated with three majorcategor- 
ies: phosphoproteins, alternative splicing, and negative 
regulation of metabolism (repressors), thus suggesting 
that the MRL/MpJ mouse is likely to produce more 
phosphoproteins, repressors, and alternatively spliced 
transcripts. Interestingly, intragenic DNA methylation 
has been associated with the use of alternative promo- 
ters and the synthesis of alternative transcripts, 40 while 
positive regulation of phosphate metabolic processes in 
human foetal heart tissues has been concluded from 
the analysis of DNA methylation microarray data. 41 

An additional term shared by the term sets obtained 
for the genes hypomethylated in liver, heart, and spleen 
of the MRL/MpJ mouse was associated with embryonic 
morphogenesis. This suggests that the DNA methylation 
pattern of MRL/MpJ mouse might have retained some 
embryonic features and the possibility will be discussed 
in the next chapter. 

The remarkablefeaturesofthefunctionalterms repre- 
senting the hypermethylated genes identified in the 
heart and spleen of the MRL/MpJ mouse were associated 
with three groups: G-protein-coupled receptors, includ- 
ing m ost ly o If acto ry rece pto rs, I g- 1 i ke V-ty pe , a nd n at u ra I 
killer receptor Ly4 9 genes. Over 1 500 olfactory receptor 
genes are one of the most abundant gene family in 
mouse, and this may explain their overrepresentation. 
Nevertheless, it is worth noting that altered olfactory 
function has been reported in the MRL mouse model 
of the central nervous system lupus. 42 The majority of 
se I ected I g- 1 i ke V-ty pe ge nes be lo n gs to t he c I uste r of se - 
lection and upkeep of intraepithelial T cell genes. The 
region of chromosome 4, where the selection and 
u pkeep of i ntraepithel ia I Tcel Is genes a re located, is poly- 
morphic, and we found that the MRL/MpJ mouse lacks 
the nucleotide sequences corresponding to the DNA 
methylation peaks in the promoter regions of Ski nt3 
and Skint4 (not shown), so it is not the epigenetic but 
the genetic impact rather which should be considered 
for this group of genes. The hypermethylated genes of 
natural killer Ly49 receptors included mostly killer cell 
lectin-like receptor subfamily A genes: Klral 2, Klral 3- 
ps, Klra4, Klral 5, Klra23, Klra22, Klral 8, and Klra33, 
located on the chromosome 6. 

None of the functional terms produced for the genes 
hypermethylated in the liver was shared by the term 
sets produced for the genes hypermethylated in the 
heart and spleen of the MRL/MpJ mouse. 

3.6. Do the genomic DNA methylation patterns of the 
MRL/MpJ mouse retain embryonic features? 
Enhanced healing typical of foetal and neonatal 
period suggests that the regeneration abilities of the 
MRL/MpJ mouse could be embryonic relics preserved 
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in the adult organism. As indicated by Naviaux et ai, 24 
the adult MRL/MpJ mouse retains a selection of embry- 
onic metabolic features. With regard to finding out as to 
whether the genomic DNA methylation profiles of the 
MRL/MpJ mouse exhibit embryonic traits, we investi- 
gated the representation of the genes associated with 
embryo development among the genes differentially 
methylated in the MRL/MpJ mouse, as well as we evalu- 
ated the overall similarity of the MRL/MpJ and embry- 
onic methylation patterns. 



3 . 7. Genes associated with embryonic development are 
overrepresented among those hypomethylated in 
the MRL/MpJ mouse 
The gene set enrichment analysis for the genes hypo- 
methylated in the MRL/MpJ showed excessive numbers 
of genes associated with embryonic development 
and homeobox genes (Table 2).'Embryonicorgan mor- 
phogenesis,' one of the terms related to embryonic 
development, was shared by the functional term 
sets produced for the three examined tissues. Three 



Table 2. Excessive numbers of homeobox and embryonic morphogenesis genes are hypomethylated in the tissues of adult MRL/MpJ mouse 



Terms identified by gene set enrichment analysis with P-value of <0.05 


Number of associated genes 


Heart 




Spleen 


Homeobox 




39* 


20* 


GO:0048562 — embryonic organ morphogenesis 


0 




9 


Developmental protein 


1 3 


70* 

79 


77* 

27 


GO:0001501 — skeletal system development 






12 


GO:0001568 — blood vessel development 


g 


30* 




GO:0001654 — eye development 




14 


n 


GO:0001655 — urogenital system development 




17 


7 


GO.0001944 — vasculature development 


s 


30* 




GO:00020S8 — lens development in camera— type eye 


3 






GO:0007423 — sensory organ development 


1 2* 


23 




GO.0007507— heart development 




Z4 




GO.0009792 — embryonic development ending in birth or egg hatching 




39 


15 


GO.0014031 — mesenchymal cell development 








GO:0021510 — spinal cord development 




9 


5 


GO.0021517 — ventral spinal cord development 




6 


4 


GO:0030879 — mammary gland development 


4 




5 


GO:0030900— forebrain development 




20 


7 


GO:0030902— hindbrain development 




9 


6 


GO:003 1 1 75 — neuron projection development 


6 




10 


GO:0035270 — endocrine system development 




12 




GO:0035295— tube development 


6 


32* 




GO:0043009 — chordate embryonic development 




39 


14 


GO:0043010 — camera-type eye development 


6 


13 


7 


GO: 00435 83— ear development 


6 




11 


GO:0048568 — embryonic orgn development 


6 


35* 


11 


GO:0048598 — embryonic morphogenesis 




45* 


14 


GO:0048666 — neuron development 


8 


13 


GO:0048706 — embryonic skeletal system development 




12 


6 


GO: 0048732— gland development 


6 


22* 




GO: 0048839— inner ear development 


5 






GO : 0060429 — epithelium development 




27 


11 


GO:0060485 — mesenchyme development 






5 


Number of genes in analysis 


180 


1063 


336 



The gene set enrichment analysis carried out using DAVID showed that homeobox genes and a number of other genes involved 
in multiple developmental processes were significantly overrepresented in the group of the MRL/MpJ hypomethylated genes 
(Supplementa ry Table S1 ).Theenrichment terms which are statistically significant with modified Fisher's exact test (p-values 
<0.05) are indicated with orange fields (darker shade of grey in the print) and those with p-values <0.05 after Benjamini cor- 
rection are indicated with red fields (darker shade of grey in the print). 
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different groups of genes associated with embryonic 
organ morphogenesis which were identified in the 
liver, heart, and spleen shared two elements: the 
insulin induced gene 2 (Insig2) and forkhead box C2 
(Foxc2) (Supplementary Table S6).The latter could be 
of particular interest of regenerative medicine as a 
potent regulator of angiogenesis. 43 

Homeobox genes have been found to be hypomethy- 
lated in the embryos of the C57BL/6J mouse in hearts 
and livers (spleens were not examined in that study). 44 
Homeobox genes were overrepresented among the 
genes differentially methylated in the spleens and 
livers, but not in the hearts of the MRL/MpJ mouse. 
Thirty-nine and 2 0 homeobox genes were found to be 
hypomethylated in the liver and spleen in the MRL/ 
MpJ mouse, respectively (Table 2 and Supplementary 
Table S6). Seven of these homeobox genes were hypo- 
methylated in both tissues. 

3.8. Genomic DNA methylation patterns of the MRL/ 
MpJ mouse reveal lower inter-tissue diversification 
To evaluate the overall similarity between the 
genomic DNA methylation profiles of embryos and 
the adult MRL/MpJ mouse, we applied Pearson's correl- 
ation test. We analysed the DNA methylation data 
obtained in our study and that reported by Liang 
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et al. 44 As we were comparing the data produced by 
different generations of Nimblegen microarrays, we 
restricted the analysis to the transcription start sites 
shared by both platforms. The analysis of correlation 
data showed that the inter-tissue diversification of the 
genomic DNA methylation profiles was lower in the 
MRL/MpJ than in the C57BL/6J mouse (Fig. 4). 
Assuming that the DNA methylation patterns reflect 
tissue differentiation, which advances with develop- 
ment, the lower diversification of inter-tissue DNA 
methylation patterns seems likely tooccuratea rlierde- 
velopmental stages. With regard to this hypothesis, we 
analysed the changes in correlation coefficients 
between the genomic DNA methylation profiles in 
embryo, newborn, and adult tissues using the data 
reported by Liang etal. 44 to find out that in most cases 
the inter-tissue correlations of DNA methylation profiles 
decrease with age (Fig. 4, inset). To our knowledge, low 
inter-tissue diversification of the embryonic DNA 
methylation profiles has never been reported to date. 

The genomic DNA methylation profiles of the 
adult MRL/MpJ mouse are not embryonic, but they 
display several features resembling those of 
embryos. The analysis of genomic DNA methylation 
patterns in the embryonic and newborn tissues of 
the MRL/MpJ mouse is necessary to provide further 
evidence. 
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Figure 4. Genome-wide DNA methylation profiles display lower inter-tissue diversification in the MRL/MpJ than intheC57BL/6J mouse. Inset: 
Similarly, we observed lower inter-tissue diversification of genome-wide DNA methylation profiles reported by Liang eto/. 44 in embryos and 
newborns than in adults (with the exception for the correlations between the heart and bra in). The inter-tissue correlations between pairs of 
genomic DNA methylation profiles were determined using Pearson's correlations coefficients. All Pearson's correlation coefficients were 
statistically significant with P-values of <0.0001. Pearson's correlation coefficients were calculated for the pairs of genomic DNA 
methylation arrays obtained in this study: heart-spleen, spleen-liver, and heart-liver for both the MRL/MpJ and C57BL/6J mouse. 
Analogically, Pearson's correlation coefficients were determined for the genomic DNA methylation profiles in different tissues of murine 
embryos, newborns, and adults reported by Liang et al. 44 The microarray data were extracted from the dataset deposited in the GEO 
database under the accession number GSE2141 5. The Pearson's correlation coefficients shown in the inset were determined for the 
average values obtained from two biological replicates, but similar trends were found in the correlations determined for individual 
animals (not shown). 
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Figure 5. Excessive numbers of high DNA methylation peaks in the 
genome-wide DNA methylation profiles in the tissues of MRL/ 
MpJ mouse. It is assumed that high methylation peaks have the 
values of 4.0 and over. 



3.9. Enhanced proportion of high methylation peaks 
in the MRL /MpJ compared with C57BL/6J mouse 
In addition to the above-discussed low inter-tissue 
diversification, another characteristic feature of the 
genomic DNA methylation profiles of the MRL/MpJ 
mouse, we found, is an enhanced proportion of high 
methylation peaks (with the values of 4.0 and over) in 
comparison with the C57BL/6J (Fig. 5). 



3.10. The balance of DNA methylation is shifted 
upstream the genes in the MRL/MpJ mouse 

It is not only the degree of CpG methylation in the 
regulatory region of a gene, but also the distance of 
methylated CpGs from the transcription start sites that 
contributes to the repression of gene expression. 39,45 In 
other words, a DNA methylation peak though mapped 
to a gene may have different impacts on gene expression 
dependent on its location. 

To examine the distribution of DNA methylation 
peaks in relation to the distance from their promoters, 
we calculated peak numbers at different distance 
ranges from transcription start sites and primary tran- 
scripts in hearts, livers, and spleens of the MRL/MpJ 
and C57BL/6J mouse (Fig. 6). The microarray design 
of the Nimblegen platform we used includes mainly 
the transcription start sites associated with the basal 
promoters. 

We found that the numberof DNA methylation peaks 
situated over 500 bp upstream from transcription start 
sites in the MRL/MpJ substantially exceeds than that in 
the C57BL/6J mouse (Fig. 6A). The number of DNA 
methylation peaks mapped within close promoter 
neighbourhood, in the range from 500 bp upstream 
up to 1 00 bp downstream from transcription start 



sites, is in turn higher in the C57BL/6J than that in the 
MRL/MpJ mouse (Fig. 6A). Such differences were 
revealed in all examined tissues and they do not 
reflect the differences in total peak numbers. The 
numbers of DNA methylation peaks located within 
the primary transcripts and up to 1000 bp down- 
stream, which are mostly the intragenic ones, were re- 
markably higher in the C57BL/6J than that in the 
MRL/MpJ mouse in the spleens and livers, while in the 
hearts they were similar (Fig. 6A). 

We also found that the numbers of DNA methylation 
peaks mapped over 5000 bp upstream and over 
1000 bp downstream from primary transcripts, we 
may consider intergenic, are lower in the MRL/MpJ com- 
pared with the C57BL/6J mouse (Fig. 6A). However, the 
microarray platform, we used, is a promoter one, and 
few probes are derived from the intergenic regions. 
Consequently, this estimate of intergenic DNA methyla- 
tion may not be representative. 

Interestingly, the differences described above are 
dramatically more pronounced for the differentially 
methylated genes (Fig. 6B). 

The results indicatethatthe balance of DNA methyla- 
tion in the MRL/MpJ mouse is shifted towards the distal 
alternative promoters and other distal regulatory ele- 
ments as, for example, enhancers. This shift is particu- 
larly remarkable for the differentially methylated 
genes. However, it should be stressed that the observed 
methylation of distal regulatory regions is not prevail- 
ing, but rather more accentuated in the MRL/MpJ 
than in the C57BL/6J mouse. Alternative promoters 
are known to be associated with development and 
tissue and cell differentiation, 46,47 which implicates 
the question as to whether this trait of DNA methylation 
profiles could impactthecellulardifferentiation poten- 
tial in the MRL/MpJ mouse. 



3.11. Genes with similar DNA methylation status in the 
heart and liver of the MRL/MpJ, but not the 
C57BL/6J, mouse 
As regeneration abilities of mammalian liver remark- 
ably exceed those of the heart, the genes that display 
strong differences in the DNA methylation status 
between heart and liver could be of particular interest 
with regard to regeneration. Therefore, we singled out 
the genes that had 'liver-like' methylation status in the 
MRL/MpJ, but not in the C57BL/6J heart. We identified 
two, relatively small, groups of such genes. One group 
consisting of 1 8 genes and CpG islands were the genes 
that were hypermethylated inthe heartand hypomethy- 
lated in the C57BL/6J liver, while hypomethylated both 
in the heart and liver of the MRL/MpJ mouse. Three 
genes involved in embryonic morphogenesis: Epha2, 
Pax2, and Gatad2a (Table 1 ) deserve particularattention 
in this group. The other set, including 47 genes and CpG 
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Figure 6. The numbers of DNA methylation peaks at different distance range from transcription start sites and primary transcripts in the MRL/ 
MpJ and C57BL/6J mouse. DNA methylation peaks were divided into four categories, according to their location from the genes they were 
mapped to: (i) distal, mapped from 5000 bp upstream up to 500 bp upstream from transcription start sites (TSSs), (ii) proximal, mapped 
from 500 bp upstream up to 1 00 bp downstream from transcription start sites, (iii) mapped within primary transcripts (PT) and 1 000 bp 
downstream from PT3' ends, except those up to 1 00 bp from TSSs, (iv) the peaks with in remote CpG islands that were not mapped to any PT, 
i.e. those located over 5000 bp upstream from TSSs and 1 000 bp downstream from PT 3' ends. (A) Peak counts including all DNA 
methylation peaks of >2.0. (B) Peak counts including differentially methylated peaks only (the peaks selected according to the criteria 
specified in the chapter: 'Genes differentially methylated in the MRL/MpJ mouse'). 
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islands, were hypermethylated in the C57BL/6J heart 
and hypomethylated in the liver, while hypermethylated 
both in the heart and liver of the MRL/MpJ mouse. Two 
groups of genes, 1 0 genes associated with acetylation, 
and 3 genes encoding LIM domains were overrepre- 
sented in this group (Supplementary Table S4). The 
result was not statistically significant for the first group, 
but it includes as an important gene as histone deacety- 
lase 1 (Hdad), which is one of key factors responsible 
for chromatin remodelling. The repression of this factor 
could prevent chromatin condensation, thus promoting 
the transcription of a number of target genes. The over- 
representation of three LIM zinc-binding domain genes 
(Csrpl, Abliml , and Fhl1) was statistically significant, 
and LIM domains are known to play important roles in 
the regulation of developmental pathways. 48 

3.12. Remarkable representatives of the differentially 
methylated genes 

As we identified as many as hundreds of differentially 
methylated genes in three tissues of the MRL/MpJ 
mouse (Supplementary Table S1), we would like to 
give a short description of a selection of these genes 
(Table 1 ). The selection includesthegenesdifferentially 
methylated in all examined tissues and several others 
that could be of exceptional interest due to their poten- 
tial importance to regenerative abilities. A few genes 
from this list which arrested our distinct attention and 
were not discussed before are discussed below. 

Since the MRL/MpJ mouse exhibitsenhanced healing 
abilities, wound repair-related genes could be expected 
to be found in the group of the genes differentially 
methylated in this strain. As the genes related to 
wound healing are not a precisely defined category, it 
is difficult to evaluate if such genes are overrepresented 
amongthosedifferentially methylated in the MRL/MpJ. 
We found two genes associated with wound healing 
among the genes hypomethylated in the MRL/MpJ 
mouse that deserve particular attention: platelet- 
derived growth factor alpha (Pdgfa) and platelet- 
derived growth factor receptor alpha (Pdgfra). The 
roles of platelet-derived growth factors (platelet- 
derived growth factor beta in particular) and their recep- 
tors in wound healing have already been established. 49 

The promoter regions of the gene encoding de novo 
methylase Dnmt3b are hypomethylated in the heart 
and spleen of the MRL/MpJ mouse. Dnmt3b could be 
of particular importance as it is expressed in embryonic 
stem cells and embryo inner cell mass and when 
expressed in adults could contribute to epigenetic 
plasticity. 

The promoter regions of the growth arrest and DNA 
damage-inducible 45 alpha (Gadd45a) gene were 
hypermethylated in the MRL/MpJ mouse in the 
spleen and heart (but not in the liver). Interestingly, 



enhanced methylation of this region has been observed 
in response to folate intervention during spinal cord 
regeneration. 2 

Regeneration abilities are often associated with the 
pluripotency of stem cells. Nevertheless, the list of 
genes differentially methylated in the MRL/MpJ 
mouse does not include the key stem cell marker 
genes, such as Oct4 (Pouf51), Klf4, Sox2, Nanog, and 
Isletl (which are covered by the array), where the 
three latter have been reported to display enhanced 
expression in the MRL/MpJ mouse. 24 

4. Concluding remarks 

The analysis of genomic DNA methylation profiles 
revealed hundreds of genes hypo- and hypermethy- 
lated in the heart, spleen, and liver of the MRL/MpJ 
mouse strain in comparison with the C57BL/6J 
control, and 37 of these genes were differentially 
methylated in the three examined tissues. The distin- 
guishing features of the genomic methylation patterns 
in the tissues of adult MRL/MpJ mouse, we found, were 
the shift of DNA methylation balance upstream the 
genes and relatively low inter-tissue diversity, resem- 
bling that we observed in previously reported 
genomic DNA methylation profiles of embryos and 
newborns. 44 The possibility that the DNA methylation 
patterns of the MRL/MpJ mouse retain embryonic 
relics is supported by another observation, which is 
the overrepresentation of the genes involved in em- 
bryogenesis among the genes hypomethylated in the 
MRL/MpJ mouse. To our knowledge, this is the first 
study to report embryonic features of the genomic 
DNA methylation profile in adult organism. 

This study presents a group of genes differentially 
methylated in the MRL/MpJ mouse, which could be 
analysed as candidate genes of regenerative phenotype. 
The examination of additional reference strains is ne- 
cessary to narrow the group by eliminating the genes 
related to inter-strain differences that are not linked 
with regeneration. Our results show that the embryo- 
like features of the DNA methylation patterns are 
likely to contribute to enhanced regenerative capacity 
observed in different tissues of the MRL/MpJ mouse. 
The analysis of the DNA methylation changes in the 
wound area in the MRL/MpJ and non-healing reference 
strains is necessary to examine the epigenetic differ- 
ences which occur in response to injury. 

The Fas ]pr mutant of the MRL/MpJ mouse has been 
widely investigated as a model of systemic lupus erythe- 
matosus and other autoimmune diseases. The MRL/ 
MpJ mouse has been also examined in the context of 
cholesterol metabolism. 38 Although we focused on 
wound healing and regeneration aspects rather, we 
think that our data could also be useful for autoimmune 
and lipid metabolism studies. 
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